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Abstract. A square lattice is introduced into the Penna model for biological aging in order to study the 
evolution of diploid sexual populations under certain conditions when one single locus in the individual's 
genome is considered as identifier of species. The simulation results show, after several generations, the 
flourishing and coexistence of two separate species in the same environment, i.e., one original species splits 
up into two on the same territory (sympatric speciation). As well, the mortalities obtained are in a good 
agreement with the Gompertz law of exponential increase of mortality with age. 
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1 Introduction 

The understanding of species formation - groups of actu- 
ally or potentially inter-breeding populations, which are 
reproductively isolated from other such groups - is still 
a fundamental problem in biology jlj. Speciation usually 
occurs when a pre-existing population is divided into two 
or more smaller populations by a geographical barrier, like 
an island, river, isolated valley, or mountain range. Once 
reproductively isolated by the barrier, the gene pools in 
the two populations can diverge due to natural selection, 
genetic drift, or gene flow, and if they sufficiently diverge, 
then the inter-breeding between the populations will not 
occur if the barrier is removed. As a result, new species 
have been formed. 

In spite of theoretical difficulties to show convincingly 
how speciation might occur without physical separation 
12] , there is an increasing evidence for the process of sym- 
patric speciation, in which the origin of two or more species 
from a single ancestral one occurs without geographical 
isolation _3^. The most straightforward scenario for sym- 
patric speciation requires disruptive selection favoring two 
substantially different phenotypes, followed by the elimi- 
nation of all intermediate phenotypes. In sexual popula- 
tions, the stumbling block preventing sympatric specia- 
tion is that mating between divergent ecotypes constantly 
scrambles gene combinations, creating organisms with in- 
termediate phenotypes. However, this mixing can be pre- 
vented if there is assortative 0] instead of random mating, 
i.e., mating of individuals that are phenotypically similar. 
It can be based on ecologically important traits such as 
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body size (as in stickle-backs) or on marker traits that 
CO- vary with ecological traits (such as coloration or breed- 
ing behavior in cichlids)jQ. 

The present paper reports on a attempt to address 
the challenging problem of sympatric speciation using the 
widespread Penna bit-string model [3IH] for age-structured 
populations, which is based on the mutation accumula- 
tion theory for biological aging. It has successfully repro- 
duced many different characteristics of living species, as 
the catastrophic senescence of pacific salmon |^, the inher- 
itance of longevity JOI and the evolutionary advantages 
of sexual reproduction jllj , as well as a particular case of 
sympatric speciation [T^. 



2 The model 

2.1 Model without lattice 

Each individual of the population is represented by a "chro- 
nological genome" , which consists of two bit-strings of 32 
bits (32 loci or positions) each, that are read in parallel. 
One string contains the genetic information inherited from 
the mother and the other, from the father. Each position 
of the bit-strings is associated to a period of the individ- 
ual's life, which means that each individual can live at 
most for 32 periods ( "years" ) . Each step of the simulation 
corresponds to reading one new position of all individu- 
als' genomes. Genetic defects are represented by bits 1. If 
an individual has two bits 1 at the i-th position of both 
bit-strings (homozygote), it will start to suffer the effects 
of a genetic disease at its i-th year of life. If the individ- 
ual is homozygous with two bits zero, no disease appears 
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in that age. If the individual is heterozygous in that po- 
sition, it will become sick only if that locus is one for 
which the harmful allele is dominant. The dominant loci 
are randomly chosen at the beginning of the simulation 
and remain fixed. If the current number of accumulated 
diseases reaches a threshold T, the individual dies. 

If a female succeeds in surviving until the minimum 
reproduction age R, it generates b offspring every year un- 
til death. The female randomly chooses a male to mate, 
the age of which must also be greater or equal to R. The 
offspring's genome is constructed from the parent's ones; 
first the strings of the mother are randomly crossed, and 
a female gamete is produced. Mm deleterious mutations 
are then randomly introduced. The same process occurs 
with the father's genome (with Mf mutations), and the 
union of the two remaining gametes form the new genome. 
This procedure is repeated for each of the b offspring. The 
sex of the baby is randomly chosen, each one with proba- 
bility 50%. Deleterious mutation means that if a bit is 
randomly chosen in the parent's genome, it is set to 1 in 
the offspring genome. However, if a bit already set to 1 is 
randomly chosen, it remains 1 in the offspring genome (no 
back mutations). 

The description given above corresponds to the orig- 
inal sexual version of the Penna model U2|j in which at 
every time step each individual of the population, inde- 
pendently of its age or current number of accumulated dis- 
eases, can be killed with a probability Vt = 1 — A^t/A^max; 
-^max is the maximum population size (the carrying capac- 
ity of the environment) and Nt is the current population 
size. This random time-dependent death, well known as 
the Verhulst factor, is introduced in order to avoid the un- 
limited growth of the population and to take into account 
the dispute for food and space. Since there seems to be 
no biological justification for considering random deaths 
in real populations, as well as a controversial importance 
of its role in the Penna model ^3], in our simulations 
we do not consider random deaths. Instead, we adopt a 
simple lattice dynamics which also avoids the exponential 
increase of the population. The details will be presented 
in the next subsection. 



2.2 Speciation model on a lattice 

In the present case each individual lives on a given site 
{i,j) of a square lattice and, at every time-step, has a 
probability p^ to move to the neighboring site that presents 
the smallest occupation, if this occupation is also smaller 
or equal to that of the current individual's site. We start 
the simulations randomly distributing one individual per 
site on a diluted square lattice. That is, if an already oc- 
cupied site is chosen for a new individual, the choice is 
disregarded and another random site is picked out. 

At any bit position a diploid individual can have n = 0, 
1 or 2 bits set. The process of sympatric speciation is now 
attempted by defining one single bit position, which we 
take as position 11, as an identifier of the species. Mating 
occurs only among individuals of the same species (same 
value of n at position 11), which means that this locus also 



defines the mating preferences. Each able female (with age 
> R) with n such bits randomly selects a neighboring able 
male with the same n value to breed. If she succeeds, she 
generates b offspring. Then she chooses at random, again 
among its four neighboring sites, a place to put each baby, 
according to the rules below. The newborn dies if it is not 
possible to find a site respecting these rules: 

1) The selected site occupation must be < 1; 

2) If the newborn has n = 0, then it can occupy an empty 
site or a site already occupied by a single individual with 
n = 2; 

3) If the newborn has n = 2, then it can occupy an empty 
site or a site already occupied by a single individual with 
n = 0; 

4) If the newborn has n = 1 , it can occupy only an empty 
site. 

Rules 2 and 3 mean that the n = and the n = 2 
populations can share the same habitat, that is, they do 
not dispute for the same food resources. Rule 4 means 
that the n — 1 population feeds at both niches, competing 
with the other two. Theses rules replace the random killing 
Verhulst factor pointed out in the previous section. 

We start our simulations only with n = individuals. 
Due to the randomness of mutations and crossover, the 
offspring does not necessarily have the same n value of the 
parents. In our model it is exactly this randomness which 
allows the emergence of new species out of the original 
one. These populations coexist in a stable equilibrium but 
without cross-mating. 



3 Simulation Results 

The simulation starts with A^o individuals, half males and 
half females, and runs for a pre-specified number of time 
steps, at the end of which averages are taken over the 
population(s). The general parameters of the simulations 
are: 

• Minimum age of reproduction R = 8; 

• Birth rate 6 = 3; 

• Mutation rate M = 1 per bit-string (or gamete); 

• Maximum number of genetic diseases T = 5; 

• Probability to walk pw = 1-0. 

Fig- n shows how the new species A''2 emerges, within 
about a hundred iterations, from the original species A'q. 
The intermediate population A^i is only about 0.50% of 
the total population. Since the rule for an individual to 
move on the lattice depends only on the existence of a 
site with an occupation smaller or equal to that of the 
current individual's site and is completely non-related to 
the individual species, the different species may bunch to- 
gether at the same site of the lattice. Then, due to the 
reproduction rules quoted above, after several generations 
we obtain a great predominance of the two non-competing 
species A^o and N2 living at the same geographic position 
(sympatric speciation). Our results with Aq = 1600 are 
confirmed by larger simulations with Aq — 100000, and 
also by larger simulations with 10^ time steps. 

It must be remarked that the assumptions we make 
concerning mating choice and conditions for a newborn to 
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Fig. 1. Time evolution of A'^o (line, original species), Ni (+, 
mixed genomes) and N2 (x, new species), for one diploid sexual 
population simulated on a L x L square lattice with L = 800 
and No = 1600 individuals. 




Fig. 2. Mortality as a function of age averaged over 20 diploid 
sexual populations simulated on a square lattice with L — 200 
and No = 800 individuals. The dashed line corresponds to the 
fit g(age) = 0.00055 exp(0. 4* age) in this semilogarithmic plot. 



survive were adopted in order to capture some features of 
field observations and laboratory experiments of species 
which seem to speciate via disruptive selection on habi- 
tat/food preferences and assortative mating. 

For instance, in a series of papers Rice and Salt [T^ 
presented experimental evidence for the possibility of sym- 
patric speciation in Drosophila melanogaster. They started 
from the premise that whenever organisms sort themselves 
into the environment first and then mate locally, individ- 
uals with the same habitat preferences will necessarily 
mate assortatively. Others examples of sympatric speci- 
ation can be found for canids ^B], lizards ^T] and pandas 
|18| . In this latter example, the Giant Panda (Ailuropoda 
melanoleuca) and the Red Panda (Ailurus fulgens) are 
vegetarian carnivores that specialize in eating bamboo in 
Sichuan Province, China. The two species share the same 
habitats and bamboo plants. Both pandas feed on the 
same species of bamboo, but specialize in eating differ- 
ent parts of the bamboo plant. The Giant Panda feeds 
more frequently on bamboo stems, while the Red Panda 
feeds more frequently on bamboo leaves [18) . In our sim- 
ulation, disruptive selection explicitly arises from compe- 
tition for a single resource (a potentially more common 
ecological situation). In this way, we may imagine, for in- 
stance, that the original population n = is vegetarian, 
and the second population n = 2 emerging out of it con- 
sists of carnivores (thus, there is no competition between 
the two different populations). However, since the indi- 
viduals with n = 1 feed from the same resources of both 
populations (n — and n — 2), this competition for food 
reduces its abundance in the system and, combined with 
assortative mating, leads to evolutionary branching. 

The situation that better fits our simulations occurs 
in the Australian Fogg Dam Nature Reserve, where data 
have been collected J.9 from three different snake species: 
water phytons {Liasis fuscus, Pythonidae) , keelbacks {Tro- 
pidonophis mairii, Colubridae) and slatey-grey snakes {Ste- 



gonotus cucullatus, Colubridae) . All are non- venous, ovip- 
arous and active foragers, but they differ considerably in 
body sizes and dietary habits. Water phytons feed al- 
most exclusively on a single species of native rodent; keel- 
backs feed primarily on frogs and slatey-grey snakes have 
extremely broad diets (reptile eggs, frogs, small mam- 
mals and lizards). According to Ref.^Hjj the population of 
slatey-grey snakes is smaller than the other two during the 
whole year. Particularly from April to May (when neither 
the rats nor the frogs are in their peaks of abundance), the 
water phyton and the keelbacks populations are almost of 
the same size while the slatey-grey snakes population size 
is around 1/7 of this value. 

To check the main responsible element for the observed 
speciation in our model, we have also performed several 
simulations considering all possible permutations of the 
four new aspects introduced into the usual aging model 
here studied: 

1) the kind of species in the beginning of the simula- 
tion (only A'o or iVi or N2); 

2) Mating preferences: The mating occurs only among 
individuals of the same species or it occurs without taking 
into account the individual's species (randomly). 

3) Back mutation: A reversal process whereby a gene 
that has undergone mutation returns to its previous state, 
i.e., if the randomly chosen position for introducing mu- 
tation at birth is the bit identifier of the species (position 
11), with a small probability pbm, a bit set to 1 in the 
parent's genome, it is set to in the offspring genome. 
However, if this bit is already set to 0, it remains in 
the offspring genome. For all the others positions in the 
genetic strand only the possibility of harmful mutation is 
considered. 

4) With or without the coexistence rule: the newborn 
with n = (n = 2) can occupy an empty site or a site 
already occupied by a single individual with n ~ 2 (n — 0), 
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or they can born if and only if there is an empty site to 
be placed. 

AdditionaUy, for aU the cases, we also analyzed if there 
would be any finite size effects related to the number of 
sites which the newborn can be placed: the newborn can 
be randomly placed only on the nearest-neighbor mother's 
site or in any site of the lattice. 

If the simulation starts only with individuals Nq or 
only with individuals A^i , and only harmful mutations are 
introduced, the speciation is observed only if the coex- 
istence rule is considered. However, for the case starting 
with only individuals N2, the time evolution of the pop- 
ulation presents only individuals N2 , since if only delete- 
rious mutations are considered one bit set to 1 does not 
change to one set to 0. That's why in this case, the spe- 
ciation could be observed only when the back mutations 
with pbm = 5 X 10~^ were allowed. Furthermore, in all 
cases starting only with species N2, the newborns suffer 
back- mutations with a small probability Pbm- 

If there is no mating preference, which means that indi- 
viduals from different species can mate with one another, 
then the speciation occurs, however the coexistence rule 
must be considered, as mentioned before. The results ob- 
tained when the newborns are allowed to occupy only the 
nearest-neighbor mother's sites and those ones when it 
could also occupy any lattice site show qualitatively the 
same behavior. It has been also noticed that the vanishing 
of intermediate individuals {Ni ) (Fig. ^ was not caused 
by the restriction that their offspring could be born only 
if there is an empty site in the nearest-neighbor moth- 
ers' site, since simulations taking into account the nearest- 
neighbor, the next-nearest-neighbor mother's site and any 
empty site of the lattice (independently of its neighbor- 
hood) did not avoid their extinction, which, in fact, was 
observed to be due to the competition with the others two 
species (iVo and iV2) for the same resources, that is im- 
posed through the prohibition of newborns A^i to born on 
sites which are already occupied. This conclusion could be 
reinforced with the findings of the following investigation: 
During a certain timestep t < ti, the population evolves 
considering the possibility to put newborns on empty sites 
and also the coexistence rule. For t > ti, when the popu- 
lation size has reached the equilibrium (constant in time), 
the coexistence rule is disregard and only one condition for 
the survival of the offspring is hold: the newborn will be 
born only if there exists any empty lattice site (not only in 
the mother's neighborhood) for it to be placed on. As we 
can see from the FigureOI when t < ti, the sympatric spe- 
ciation occurs and the results are similar to one showed in 
Figure^ (^0 — ^2) >> ^i- However, for t >ti, since the 
population size of the mixed species is larger than those 
with n — 2 {Nq > Ni > iV2), it is not possible anymore to 
affirm the existence of speciation, even the coexistence of 
two separate species in the same environment (site). Based 
on this result, we can conclude that the sympatric speci- 
ation obtained in our simulations is essentially caused by 
non-existence of competition between the newborns with 
A^'o and N2 and not by the finite size effects related to the 
neighborhood considered to place the newborns A'^i. 




Time 



Fig. 3. Time evolution of A^o (line, original species), Ai {+, 
mixed genomes) and A2 (x, new species) , for one diploid sexual 
population simulated on a L x L square lattice with L — 800 
and Ao = 1600 individuals, when two different strategies for 
the placement of newborns are considered. For t < ti, new- 
borns can be placed on empty sites and those ones with n — 
and n = 2 can share the same lattice site, when t > t\ new- 
borns can born if and only if there is any empty site on the 
lattice, fi = 6000. 



If the females mate with males from any species ran- 
domly, the coexistence rule is the only condition required 
for ocurring speciation, excepting when the simulation 
starts with species A''2 we must also take into account 
the back-mutations, for the reason previously mentioned. 
From this result, we can conclude that the mating pref- 
erence condition is not necessary to be imposed in or- 
der to obtain sympatric speciation in our model, but the 
species's ecological behaviour introduced through the co- 
existence rule. This conclusion is more reinforced with the 
results from all the simulations discussed before in which 
the coexistence rule always must be considered to lead the 
population to speciation. 

As a final study, we examine our populations mortal- 
ities. In 1825, based on observed death and population 
records of people in England, Sweden, and France between 
ages 20 and 60 in the nineteenth century, the British ac- 
tuary Benjamin Gompertz derived a simple formula de- 
scribing the exponential increase in death rates between 
sexual maturity and extreme old ages |20| . This formula, 
(7(age) = A X exp(b x age), is commonly referred to as 
the Gompertz's law of mortality. As Fig. [Slshows, our re- 
sults for the mortality above the minimum reproduction 
age i? = 8, are in a good agreement with the Gompertz 
law. 



4 Conclusions 

In conclusion, the results presented here are based on a 
very simple assumption that a single locus in the indi- 
vidual's genome determines the ecological behavior of the 
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individual and identifies its species, in the framework of 
the Penna bit-string model. Despite its simplicity and lim- 
ited applicability, our results clearly show the emergence 
of sympatric speciation in diploid sexual age-structured 
populations of individuals that are distributed on a square 
lattice. The introduction into the model of an specific gene 
responsible for the ecological behaviour of the individual: 
individuals with n = I and n = 2 do not compete for food 
resources, so they can bunch together at the same lat- 
tice site, tailored to reproduce recent observations of ex- 
isting species 21| that have led to the suggestion that eco- 
logical adaptation is the driving force behind divergence 
of populations leading to speciation: a gene, desaturase2, 
of Drosophila melanogaster, which confers resistence to 
cold as well as susceptibility to starvation, underlies a 
pheromonal difference and contributes to the reproductive 
isolation between some Drosophila species: the Zimbabwe 
and Cosmopolitan races. 

Since the speciation discussed in this paper is triggered 
by interaction or competion between organisms, and not 
merely by mutation, the process is not so much random 
as deterministic. In fact, the speciation process occurs ir- 
respectively of the adopted random number in the simu- 
lation, as well without assuming mating preference. Ac- 
cording to our scenario, the coexistence of the two groups 
(A'o and N2) is the main mechanism that often leads the 
population to speciation. It should be also noted that the 
change in genotypes occur within few generations. The 
speed of genetic change, of course, depends on the mu- 
tation rate, but the present mechanism is found to work 
even for any smaller mutation rate (say T < 5). 

The relationship between our model and previous phys- 
ics models, like directed percolation or reaction-diffusion 
problems, as suggested by an anonymous referee, unfortu- 
nately remains to be elucidated, even though some rela- 
tions to reaction-diffusion systems can be e stablished, a 
priori. By the way, further research is needed to explore 
that question, which it is not the goal and the context of 
the present paper. 
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